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Abstract 


Thermocapillary stability characteristics of a horizontal liquid layer heated from below 
rotating about a vertical axis and subjected to a uniform vertical magnetic field are ana- 
lyzed under a variety of thermal and electroraagnetic boundary conditions. Results based on 
analytical solutions to the pertinent eigenvalue problems are discussed in the light of 
earlier work on special cases of the more general problem considered here to show in parti- 
cular the effects of the heat transfer, nonzero curvature and gravity waves at the two-fluid 
inferface. Although the expected stabilizing action of the Coriolis and Lorentz force fields 
in this configuration are in evidence the optimal choice of an appropriate range for the 
relevant parameters ;s shown to be critic-lly dependent on the interfacial effects mentioned 
above . 


Introduction 


In recent years there has been a resurgence of interest in understanding the origins and 
possible means of controlling convective instability, especially in configurations relevant to 
material sciences in general and material processing in particular within the framework of 
the current space programs. In this context some of the basic aspects of this problem area 
have been under investigation^"'* by the present author. The contribution to be presented here 
is part of a continuing effort at the DFVLR to analyze some of the basic fluid dynamic 
aspects relevant to the material science configurations, especially in the context of space 
experiraents under reduced gravity conditions and the related ground based research. 

Since references give the general background and motivation for the particular problem 
considered here and cite the relevant literature, we shall restrict ourselves here only to 
a resport of some of the recent results obtained and discuss them in the light of those avail- 
able in the literature. While references 1-4 deal exclusively with the zero gravity situa- 
tion, we consider here specifically the simultaneous action of surface tension and gravity 
in this classical Benard - Marangoni configuration. 

Formulation of the problem 

We consider an infinite, horizontal, Boussinesq liquid layer of mean thickness d rotating 
about a vertical axis at a constant angular speed f) and subjected to a uniform magnetic induc- 
tion field of strength Bq under various typical boundary conditions to be detailed later. 
Figure 1 illustrates the configuration schematically and is followed by a list of the symbols 
for dimensional quantities occurring in the later development. The details of the formulation 
incorporate the features introduced by Scriven and Sternling ^ and Smith ° extending the 
pioneering work of Pearson’. 
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Figure 1. The Bfenard - Marangoni configuration 


List of symbols 


B = Coefficient of thermal 1 3p 

volume expansion, ^ 3T 

Y = Electrical conductivity 

fiT = Applied temperature difference (To -Ti ) 
(g * Amplitude of the disturbance wave at 
the two-fluid interface 
= Magnetic diffusivity (YVm)”' 

K = Thermal diffusivity, K/pCg 

X = Disturbance wavelength, 

y = Dynamic viscosity 

Wm = Magnetic permeability 

V = Kinematic viscosity 

p = Density 

0 = Interfacial energy at the two-fluid 

interface 

n = Angular speed of rotation 


S ■ Magnetic Induction field 

Bg > Magnitude of the applied B-fleld 

K » Thermal conductivity 

T “ Temperature 

Cg - specific heat 

d > Mean thickness of the liquid layer 
g = Acceleration due to gravity 
h = heat transfer coefficient at the 
disturbed interface 

kx,ky • Disturbance wave numbers in the 
x-y-directions 

p = Time constant in the exponential 

growth/decay factor of a disturbance 
normal mode 


The liquid layer has nominally constant temperatures Tg , T^ (Tg >Tj ) respectively at its 
lower ’•■'d upper horizontal boundaries. For the sake of definiteness and simplicity the 
characteristics of the adjoining media are somewhat idealized. They are specified for the 
three cases (J) , (2) > @ of the boundary conditions (b.c.) as follows. 


In b.c. (T) we take the bottom boundary as a thermally and electrically perfect solid 
conductor. In b.c. @ the bottom boundary is a thermally and electrically perfect Insulator. 
In both cases the upper adjoining medium is taken as an electrically Insulating gas extending 
in the z-direction to infinity. The heat transfer to the gas from the liquid layer can be 
simplified (without going into the details of the possible flow in the upper medium) in terms 
of an effective heat transfer coefficient h(T)^'^ for the two-fluid interface. A detailed 
discussion of this simplification was given by Pearson’. In b.c. (3) we consider the situa- 
tion where the same ambient gas is present cn both sides of the liquid layer. 


The onset of convective instability in such a liquid layer with an initially uniform 
linear temperature profile can be formulated as a linear eigenvalue problem for the distur- 
bance amplitudes of the flow variables using the standard normal modes procedure'. We non- 
dlmensionalize the problem using ' d, d’/v, ic/d, ic/d’ , 4tYxBg/d, aT respectively as the reference 
quantities for length, time, velocity , vorticity, electric current density and temperature. 

The stability of the configuration with respect to an inf initasimal normal mode of distur- 
bance may then be stated in terms of the following eigenvalue problems in dimensionless form. 


(D’-a’ ) (D'-a’-pj) (D’-a’-pi )W - Ta(D’-a’-p 2 )DZ - Q(D’-a’)D’w = Ra • a’ (D’-a’-pj ) 0 (1 ) 

(D’-a’-pj) 0 + W = 0 (2) 

(D’-a’-pi ) Z + Ta • DW + QDX = 0 (3) 

(O’-a’-pj)X + DZ = 0 (4) 


where D ; (1/d) • d/dz and W, Z, X, 0 are respectively the dimensionless disturbance ampli- 
tudes of the z-components of velocity, vorticity and electric current density and of tempe- 
rature. 


The boundary conditions are to distinguish not only between cases (T) , @ , (5) specified 
earlier but also as to whether the neutrally stable oscillatory (p^ ^ o) or stationary (pi* 0) 
modes are considered while determining the stability boundary for the configuration. 


(a) Neutral modes oscillatory (pj /O) 


B.c. 0 

W(0) 

B.c. 0 

W(0) 

B.C. 0 

W(0) 

For Nu . ^ 

» 0 


0 = DW(0) = 0(0) = Z(0) = DX(0) 

0 = DW(0) * D0(O) = Z(0) = X(0) 

PjCg (Kinematic condition at the two-fluid interface) 


Cr 

(Bo+a^ ) 


(D’-pi-3a’)DW(0) } - i " * 9(0) } » 0 


(5) 

( 6 ) 

(7) 

( 8 ) 
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(9) 


(D’+aMW(0) 


■Ma De(0) 
NU 


For Nu = 0 
D6(0) = 0 


( 10 ) 

( 11 ) 

( 12 ) 


i^ § 57 iT 7- ) HD'-3a’-Pi)DW(0) } - W(0) . 0 

Pj (D’ +aMw(0) ♦ Ma • a- {p 3 0(O)-W(O) |=0 
B.c. at 2 = d for cases 0 , 0 , 0 are of ttie same form as thiose for case 0 at z = 0. 

(b) Neutral modes stationary (pj = 0) 

T)ie conditions (11), (12) r. jove are to be replaced by 

W(0) = 0 (which covers also (7) above) (13) 


D^W(O) + a^ • Mae (0) + { -D’W(O) ♦ 3a- DW(0) } = 0 


(14) 


Again the b.c. at z = d are of the same form for cases 0, 0, 0 as those for case 0 

at 2=0. 

The dimensionless numbers occurring in the above formulation are bo = pad' /a (Bond) , 

Cr = Kp/od (crispation) , Ma = | ( 3o/3T) AT' /yic (Mara- :oni) , Nu=dh/KAT (Nusselt) , Pr = v/ic 
(Prandtl) , Pr^ = v/nm (magnetic Prandtl) , Q=Bjd^y/ii (Chandrase)char) , Ra = ggATd’ /vk (Rayleigh) 
Ta = 2i1d’/v (Taylor), a= 2nd/i (disturbance wave number), p, =pd^/v (frequency factor for 
oscillatory disturbance mode), Pj = Pr^i'p^ , Pj = Pr*pj . The last four parameters are characte- 
ristics of a disturbance normal mode in the hydromagnetic therraocapi llary stability discus- 
sion of the configuration whereas t^e first nine describe the basic configuration. 

Briefly''* the boundary conditions (5) , (6) state the no-slip condition and the thermal and 
electromagnetic properties associated with the boundaries whereas (8), (9) cover the require- 
ments*'* of stress balance along and normal to the two-fluid interface incorporating the 
thermocapillary terms and also talcing into account the nonzero i ntorfacial curvature (Cr) , 
gravity waves (Bo) and heat transfer contribution at the disturbed interface (Nu) . 

The existence of oscillatory modes in this configuration especially at large Ta is well- 
)cnown for the buoyancy-driven case * and •■•as also demonstrated in the surface tension-driven 
case. Tlie oscillatory modes become important at low Pr but it was found ' that at least for 
Bo = 0 the incipient instability is stationary rather than oscillatory since the correspond- 
ing critical Marangoni n^imber is higher than that for the stationary mode which is independ- 
ent of Pr. It turns out that for small Bo / 0 the critical Ma^ tends to decrease and a^. -► 0 
with large Ta whereas the oscillatory modes were shown by asymptotic analysis * to occur at 
large Ma^ Ta >■' 1 , as a short wave instability with av . Thus we have some plausible 
evidence to suppose that in this configuration, where the effects of the magnetic field 
(Q(<0) which inhibits the onset of buoyancy-driven oscillatory modes (for Pr > Pr„, ) ' are 
also included, the stationary modes precede the oscillatory ones at onset of instability 

Since the practical interest in the present inves..igaticn lie.s ultimately in the suppres - 
sion of convective instability ^ consider here the case p^ = 0 in the following. It, 

however, the solution of the complete eigenvalue problem with Pj ^ 0 posed above does lead 
to oscillatory modes we have then only to compare the corresponding minimum critical 
Marangoni number with Ma^ computed here. .Ma is in any case an upper bound for stability of 
the configuration. 

The stationary modes of convective instability are given by nontrivial solutions to the 
homogeneous boundary value problems given by (1) - (10), (13), (14) for the different cases 
0 ' 0 , 0 . The secular conditions for the existence of nontrivial solutions to the 

respective homogeneous boundary value problems have been obtained by using the exact analy- 
tical solutions (combinations of trigonometric and hyperbolic functions) of (1) - (4) in the 
appropriate boundary conditions. The neutral stability characteristics of the configuration 
arc then analyzed from the resulting transcendental secular relationship in terms of the 
dimensionless parameters of the problem. Tince we have a large number of dimensionless groups 
here, We shall have to choose a suitable range of their values with some class of applications 
in view. As indicated in references *“'■ the interface curvature effects are already in evi- 
dence for such small values of Cr = 10”’, 10"'' yielding stability characteristics quite differ- 
ent from those for Cr=0. Using the thermophysical property data available in the literature 
the parameter ratio Bo/Cr = gd’/vic at g= 9.81 m/s' with d = mm has values of 0(10' ) as 
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shown In Table 1 for some substances of interest 


f 


i 


i 


Silicone oil 
(Dow-Corning 
200) 

Cu-melt Al-Cu-melt 

Ga As-melt 

Si-melt 

4.5* 10* 

6.07'* '10* 5. 1 * 10* 

7.1*10^ 

1 .82 * 10* 


The corresponding values for different levels of gravity and the oprlate size of the 
layer for experiments In space missions can be estimated from Table 1 . We can also use them 
for estimating parameters su-'h as Ra and Bo say by choosing Cr ' 10''* , 10“" to demonstrate 
the effects of nonzero interracial curvature. It is found that Cr « 10”’ , 10“'' for Silicone 
oil (Dow-Corning 200) , used frequently for convection experiments, when d«mm, cm respec- 

tively. Thus for experiments in a terrestrial laboratory in the mm size and in an orbital 
laboratory in the cm range can be covered by considering Bo « 0.05, 0.5 and Cr ■ 10”* , 10"'' 
to emphasize the effects of interfacial waves. Table 2 gives some representative values for 
Q and Ta. 


Table 2. 

Values of Q 

at Bg >0.5 tesla and Ta at 0 

= 500 rpm for d “ mm 


Silicone oil 
DC 200 

Ai-meit Gj As- si-melt 

melt 

Cu-melt 

C 


9.02 * 10* 5.57 * 10* 1 .14 * 10* 

8.43 * 10* 

lA- 



1.95 * 10* 3.37 * 10* 2.98 * 10* 



Results and discuss ion 


The eigenvalue relationships giving the stability characteristics of the present configu- 
ration under b.c. (?) , 0, 0 have been obtained by investigating various special cases: 

Ta = 0, Q ^ 0; '' Ta * 0, Q ■ 0 : * Ta,Q » 1 ; ^ Ta , Q = 0, pi * 0 * all with P.a « 0 • Bo i.e. , under 
zero gravity bringing out the essential differences between Cr = 0 and Cr / 0 . 




Ml MAmu •* rtlilve A • I Ta • M 






Figure 3. Variation of Mac with rotation 
(Ta) and magnetic field (Q) neglecting 
Interfaclal effects (Cr = 0 • Nu) under zero 
gravity for b.c. (T) , (|) , @. 


Figure 2. Neutral stability curves for the onset 
of thermocapillary convective instability in a 
liquid layer under zero gravity; Effects of inter- 
facial curvature (Cr) and heat transfer (Nu) un<*.er 
the influence of (a) rotatifln alone and (b) mag- 
netic field alone for b.c. 0 . 

The unstable domain is above the respective curves. 
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Figure 2 shows the neutral stability curves of the present configuration under the action 
of (a) rotation alone and (b) magnetic field alone for different Nu, Cr at Bo < 0 - Ra for b.c. 
(T) . We notice first of all the radical departure of the stability characteristics for Cr K 0 
from those of Cr > 0, namely, that there exists strictly speaking no minimum critical Maran- 
goni number when Cr / 0 as was first shown by Scriven and Sternling ’ for Ta = 0 » Q. Asymp- 
totic analysis in the limit a -► 0 shows that Ma f (Ta,Q) (Nu+ 1 ) /a’ ,a -* 0 for Cr ■= 0 whereas 
Ma g (Ta,Q) (Nu+1 ) a’ /Cr , a-^0 for Cr / 0 under b.c. (l) . The numericai results shown confirm 
this limiting behaviour as well (note the linearity of the curves for small a) . The above 
formulas Incidentally include the factor (Mu-fl) missing In those of reference S (Table 1, 
p. 333) for Ta = 0 = Q. 

For sufficiently small Cr/g(Ta,Q) we may speak of r quasi-critical Marangoni number Mag 
which is approximately equal to that calculated using Cr ■ 0 in earlier literature**'*’ . Since 
the unstable long wave band increases in size with Ta and Q (for Ta >vl , Q << 1 j Ta«1, Q »1 
respectively the corresponding band widths are 0 (Cr 'Ta) and 0(Cr/7), Cr must Indeed accord- 
ingly be smaller for this approximation to hold at higher Ta , Q . As shown In Figure 2 the 
effects of heat transfer (Nu/0) at the two-fluid Interface are stabilizing in that the 
unstable domain Is pushed upward along the Ma-axls with Increasing Nu. 

Figure 3 shows the monotonically Increasing stabilization pot.'ntially to be achieved by 
Increasing rotation (Ta) and magnetic field (Q) under the three typical b.c. (T) , (2) , Q) 
which. It may be noted, are In decreasing order of stability amongst themselves. The results 
shown agree with those in references 12,14 for Cr= 0. Asymptotically Ma^ ’ 0(Q) for Q» 1 , 

Ta < < 1 and Ma^, =0(Ta) for Ta » 1 , Q « 1. Note that the asymptotic range Is attained faster 

by Ta than by Q due to the influence of rotation on the flow field in general and vorticity 
in particular. The differences between the b.c. per.nlst longer In the case of magnetic field. 
The situation is analogous In the case of buoyancy * 

Apart from their formal interest the results shown in Figures 2 , for Cr = 0 may also be 

seen as useful approximations for sufficiently i-mall Cr and at low levels of gravity provided 

the long wave instabilities are considered relatively harmless. The relevant ranges of the 
parameters will become apparent in the later discussion. 




Figure 4. Neutral stabi.l- Figure 5. Correlation of the critical Marangoni 

ity curves for b.c. (T) with number Ma^. with Bo/Cr for b.c. (J) with Nu = 0, 

Nu = 0. Bo * 0. 05 at low grav- Ra = 0. 1 at Cr « 10~* , 10”" ; Bo = 0 , 01 , 0.05 

ity (Ra*0.1) for different 
Ta (Q * 0. 1 ) and Q (Ta= 0. 1 ) 


Figure 4 shows the further departure of the stability characteristics of the configuration 
from those at Cr * 0 when we consider low Bond number and crlspatlun effects together. We no- 
tice first of all the reinstatement of an absolute minimum ciitlcal Mac onset of 

convective instability as was first shown by Smith* for Ta » 0 ■ Q . The neutral stability 
curves for Cr • 10"*show two minima, one at a « 0 and the othvr at finite a. Even for Cr •- 10"* 
the same feature can be reproduced at Bo >0.5. This Is due to the fact that the long wav; 
stability characteristics depend on the ratio Bo/Cr and not individually on Bo,Cr. The 
occurrence of double minima has been confirmed for Bo/Cr* 200. The lesser of the two minima 
is then the critical Mac onset of instability. For small Ta(<25.8) and small Q(<18) 

we observe that the critical wavelength 1^. corresponding to Ma_ is finite whereas at higher 
Ta and Q, l^c Infinite at onset of instability. It may also be mentioned that when is 
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finite it does correspond to the value for Cr= 0. Ma^. clearly decipeases roonotonically with 
larger Ta,Q and the corresponding is then infinite. Thus, allowing for gravitational 
waves and crispation effects leads to long wave instability at a low but finite Ma^ for 
large Ta and Q. 

Figure 5 shows the correlation of Ma^ with Bo/Cr for b.c. (l) with Nu= 0, Ra= 0.1 at 
Cr = 10"’, 10"“* ; Bo = 0.01, 0.05. Along the continuous parts of the curves Xj. is infinite and 
along the broken ones X^ is finite. The lattjr situation is found to occur at low Ta(<25.8) 
and Q < Q* (Q* = 18, 12, 8.5 respectively for Ta = 0.1, 15, 20) and large enough Bo/Cr. The last 
provision is to be recognize I along the curves for Ta < 26 where the respective curves split 
off at Q* into two branches applying separately for Bo =0.01 (Xj,-^~) and Bo =0.05 (X(. finite) 
eventhough both correspond to the some value of Cr = lO"’. For large enough Tr and Q values 
the correlation with Bo/Cr is universal and a(.= 0. Furthermore we notice that at large 
Q(''>800) all the curves for Ta £ 500 merge. This implies a certain "saturation effect" as far 
as the influence of rotation is concerned while acting together with the magnetic field as 
a stabilizing agent. 





Figure 6. Variation of Ma_ with ro- 
tation (Ta) and magnetic field (Q) 
for b.c. (T) with Nu = 0 

Figure 6 shows the variation of Mac with Ta and Q (225) wherein the monotonic decrease of 
Ma^. is to be noted even for Bo/Cr = 500 in contrast to the initial increase observed in Fig- 
ures 4, 5 for lower Ta, Q. (In Figures 6-3 the respective constant parameter values are 
indicated in the inset.). The computations snow that Maj, hardly changes with Ra. In fact 
for Cr = 10”’, Bo = 0.05 the upper right quadrant of Figure 6 shows that even up to Pa = 1000, 
Ma^ is that given by the long wave limit (Malg-Q . However for higher Bo/Cr (=500 shown) the 
buoyancy effects become noticeabl" from Ra > 300 for Ta =0.1, 10 at Q= 25 since now Q* in- 
creases for Ta= 0.1, 10 from 18, 12 respectively at low Ra (=0.1) to 47, 46 at moderate 
Ra (=1000). This is indicated in the lower right hand quadrant of Figure 6. This latter 
range, where buoyancy effects become noticeable, is distinguished by the broken curve along 
which Xj, is finite. (The curves for Ta=0.1, 10 are hardly to distinguish on the scale drawn 
but they end, when extended, respectively at approximately Q, = 16.5 and 13.2 on the Q-axis.) 

Figure 7 shows the correlation of Ma^. with (Bo/Cr) (Nu+1) f (Ta, Q) for different combina- 
tions of Bo, Cr, Nu. The coefficient functions f (Ta, Q) shown have been confirmed numeri- 
cally for various combinations of the pareuneters as long as the buoyancy effects are not 
noticeable. The "universality" of these correlation functions depends slightly on the para- 
meter range but is found to be within a few percent at a = 0.02 chosen to represent the limit 
a 0 . Another feature to be noted from Figures 4 - 7 is that a^ -»■ 0 as Ta, Q increase and 
3c = 0 for all Ta, Q greater than some not too large a value. This is in contrast to the 
common finding of the earlier studies (wherein Cr was set equal to zero a priori), 
namely, that a^, increases with Ta and Q. Here we see that as long as the buoyancy effects 
do not dominate, the stationary form of instability sets in only at 3^=0 for sufficiently 
large Ta and Q. 

Now we turn to the effect of the boundary conditions on t1^ stability characteristics of 
the configuration. In all the three cases of b.c. (J) , (2), @ the seune trends in the varia- 
tion of Mac observed for low Ra. Ma^, is proportional to Bo/Cr and decreases monotonlcally 


Figure 7. Variation of Mag with Bo, 
Cr, Nu for b.c. (T) 
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with Ta, Q for low Bo/Cr as demonstrated in Figure 8(a) for Bo/Cr = 10, 100. 



Figure 8(a). Effect of boundary conditions on thevariation of Ma^. at Bo/Cr = 10, 100 



Figure 8(b). Effect of boundary conditions oti the variation of Ma^ at Bo/Cr = 500 


We also note that at low Ta and 0, Ma^, for b.c. (l) is higher that that for b.c. Q) and the 
latter in turn is higher than that for b.c. @ . This indicates the decreasing degree of 
stabili^ Imparted by the degrees of freedom allowed by the three types of boundary condi- 
tions (5 ' © ' (|) that order. This feature is similar to that for the buoyancy-drive.' 

convective instability* although the laoundary conditions there are different. For large Ta 
and Q, however we notice ( cf. Q >500, or Ta = 500) the boundary conditions can no longer be 
distinguished from each other. In the present case the role of the b.c. is further enhanced 
via the dependence on Bo/Cr. The lower Bo/Cr, the lesser is the influence of b.c. even at 
low Ta,Q (cf. Ta = 1,50 for Bo =0.01, Cr = 10~* shown by dashed curves in the lower part of 
Figure 8 (a) . ) 

Complementary to the results in Figure 8(a) those in 8(b) demonstrate that for larger 
Bo/Cr ( = 500 for three different combinations of Bo,Cr) a more pronounced effect of the 
boundary conditions on the variation of Mag and in particular that Ma^, can decrease as well 
as increas'* with Ta and Q depending on the range of parameters. Again at large Ta, Q the 
distinction between the boundary conditions decreases. 

Conclusions 

The onset of stationary convective instability driven by both density- and surface-tension- 
gradierts in a horizontal liquid layer heated from be ow can be suppressed by means of rota- 
tion about a transverse axis and by a transverse magnetic field. But the stabilizing influ- 
ence of these two agencies is subject to considerable qualifications in view of the effects 
of curvature and gravity waves at the two-fluid interface. The larger the ratio Bo/Cr, the 
greater the range of stabilizing action in terms of Ta, Q for all the boundary conditions 
considered and relatively greater for b.c. Q) and (5) than for b.c. ©. The Influence of 
the individual b.c. (T) , @ , (3) becomes in3istinguTshable at larger Ta, Q and at lower 
Bo/Cr. Since Mac decreases with Ta, Q (for sufficiently large Ta, Q) and ac -» 0, an optimal 
parameter range for the combined stabilizing action of rotation and magnetic field must be 
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sounltt appropriately. Allowance for heat transfer he ambient gas is generally stabilizing. 

In the low gravity situation (small Bo and Ra < 100) the buoyancy effects do not percepti- 
bly influence t)ie onset of instability except at low Ta and Q. In this range the onset of 
instability is at a finite wave number a^ / 0 whi^'h is independent of Bo, Cr as may be expect- 
ed. (Ma^. corresponds otherwise to = 0 and as shown « Bo/Cr for Nu = 0.) The general 
problem of interaction between buoyancy and surface tension will be considered in a later 
report but it seems legitimate to draw a partial conclusion on the basis of results shown 
here, namely, that the buoyancy-dominated situation tends to prefer finite wave length in- 
stability while the capillarity-dominated situation including the effects of interfacial 
curvature and interfacial gravity waves tends to favour the infinitely long wave mode of 
instability. This conclusion is qualitatively in constrast to that of earlier studies on this 
configuration ignoring the interfacial effects altogether (Bo = 0 = Cr) . This stems only from 
the nonzero Bo/Cr and does not explicitly depend on the (finite) value of the mean surface 
tension 

The question of oscillatory modes of instability has been bypassed here on the basis of 
asymptotic results indicating that the incipient instability is stationary for large Ta(Q=0). 
The results for the finite range of Ta and Q need of course to be examined in order to 
confirm whether Ma^. calculated here is indeed the absolute minimum critical Marangoni number 
for the onset of instability. 
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